From ballistic to Brownian vortex motion in complex oscillatory media 
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We show that the breaking of the rotation symmetry of spiral waves in two-dimensional complex 
(period-doubled or chaotic) oscillatory media by synchronization defect lines (SDL) is accompanied 
by an intrinsic drift of the pattern. Single vortex motion changes from ballistic flights at a well- 
defined angle from the SDL to Brownian-like diffusion when the turbulent character of the medium 
increases. It gives rise, in non-turbulent multi-spiral regimes, to a novel "vortex liquid" . 

PACS numbers: 82.40.Ck,45.70.Qj,89.75.Kd,82.40.Bj 



Chemical waves in reaction-diffusion systems with lo- 
cal excitable or simple oscillatory dynamics have been in- 
vestigated extensively both theoretically and experimen- 
tally because of their relevance for a variety of physi- 
cal, chemical, and biological processes III In (quasi) 
two-dimensional situations, spiral wave patterns are es- 
pecially prevalent and important. They determine the 
characteristics of processes such as surface catalytic oxi- 
dation reactions 3], contraction of the heart muscle 
and various signalling mechanisms in biological systems 
to name only a few examples. The dynamics of the 
core or "vortex" of spiral waves plays an important role 
in many of these phenomena: some mechanisms for spi- 
ral breakup arise from core motion (see, for example, 
Ref. 0) and moving spirals have been suggested to be 
responsible for some cardiac arrhythmias 7]. However, 
whereas the meandering instability of spiral cores in ex- 
citable media is well known 8], studies of vortex motion 
are still few in the oscillatory case, being limited to the 
unbounded acceleration characteristic of the core insta- 
bility 1^ or to erratic motion induced by spatiotemporal 
chaos or external noise 0, • 

When the local oscillations are not simple but possess 
a complex periodic or chaotic character as observed, for 
example, in chemically reacting systems 12K spiral waves 
contain synchronization defect lines (SDL) separat- 
ing domains of different oscillation phases and across 
which the phase changes by multiples of 27r. SDL have 
been observed in experiments on Belousov-Zhabotinsky 
reaction. 14, While their origin and classification 
have been investigated, little is known about how they 
influence the spiral core to which they are connected. 

In this Letter, we show that in complex oscillatory me- 
dia the emergence of SDL is accompanied by spiral core 
motion. The SDL breaks the rotational symmetry of the 
spiral wave giving rise to a generic mechanism for core 
motion that differs from the instabilities that cause the 



non-saturating core instability in simple oscillatory me- 
dia and the meandering instability in simple excitable 
media. We show that spirals in complex oscillatory me- 
dia move ballistically in directions that are determined 
by the SDL attached to the cores. In the regime of 
SDL-mediated turbulence where SDL are spontaneously 
generated, the core motion is more complicated leading 
to randomly-oriented flights of random duration. With 
stronger turbulence this dynamics leads to vortex Brown- 
ian motion characterized by a well-defined diffusion con- 
stant. We finally show that multi-spiral configurations in 
the locally non-turbulent regime lead to a spatiotempo- 
rally chaotic "vortex liquid" . 

Consider the reaction-diffusion system 

dtciv, t) = R(c(r, t)) + W'c(r, t), 

where c(r, t) is a vector of time-dependent concen- 
trations at point r in a two-dimensional domain of 
length/diameter L 16], D is the diffusion coefficient 
(taken to be the same for all species) and where the 
local kinetics is specified by the vector-valued function 
R(c(r, t)). As a paradigmatic example of a system with 
complex local dynamics we take R(c) to be given by the 
Rossler model ^3 ^i^^ = —Cy — Ry = Cx -\- 0.2cy, 
Rz = CxCz - Ccz + 0.2. For C G [2.0,6.0], the medium 
undergoes period-doubling bifurcations transforming the 
local dynamics from simple oscillatory to period-doubled 
orbits to chaotic dynamics. The rotational symmetry of 
spiral wave patterns is then broken by the ineluctable 
appearance of SDL 13]: In a simple oscillatory medium, 
one turn of the spiral wave corresponds to one period 
of the local oscillation in phase space (change of 27r in 
the phase). When this spiral wave pattern undergoes a 
period doubling bifurcation, the period of the local or- 
bits doubles (change of An in the phase). Continuity in 
the medium then forces the emergence of a narrow strip 
connected to the core, the SDL, across which the phase 





Figure 1: (color online) Left: Snapshot of the scalar field 
Ac4r,t) = l/r /q" \c4r,t + t') - c^r.t + 1' + T)\dt' in the P2 
regime for C = 3.5, L = 256, D = 0.4 (kept fixed through- 
out). The period of the PI oscillation is r = 5.95. Black 
corresponds to Acz ~ and indicates a SDL originating from 
the spiral core. Right: Magnification of the rectangular region 
in the left panel showing the SDL (red lines) at two times, a 
is defined as the angle between the spiral core's trajectory 
(thick black line) and the attached SDL. 
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Figure 2: (color online) Dependence of the angle a (filled 
diamonds) and the core velocity, v (circles) on the Rossler 
parameter C. The solid lines are guides to the eye. From left 
to right, the vertical lines mark the first period-doubling tran- 
sition to the complex oscillatory regime and the transitions to 
the turbulent regimes (see text for details). 



jumps by 27r. Two oscillations are needed to return the 
medium to its original state. Along the SDL, the local 
dynamics is period- 1 (PI) while in the remainder of the 
medium it is period-2 (P2). An example of a PI SDL 
in a P2 medium is shown in the left panel of Fig. ^ 
the period-4 (P4) regime, two types of SDL are possible: 
lines where the phase jumps by 27r (local PI dynamics) 
and lines where the phase jumps by Air (local P2 dynam- 
ics). More generally, in a medium with period 2^ dynam- 
ics 2"^"^ SDL with periodicities 2^, {k < n), may exist, 
although for the Rossler medium only the P4 regime is 
observed before the local dynamics becomes chaotic. In 
these chaotic regimes the PI and P2 SDL persist. Deeper 
in the chaotic regime turbulent states are found where 
SDL are spontaneously created and annihilated. 

First, we consider a single spiral in a disc-shaped do- 
main with no-flux boundary conditions and focus on the 
complex oscillatory regime where no PI SDL are spon- 
taneously generated (3.03 < C < 4.557). Our extensive 
simulations revealed that the simplest possible of the al- 
lowed configurations, where a single PI SDL is attached 
to the core, is the unique asymptotic solution, even in 
the P4 and chaotic regimes where P2 SDL exist. The 
presence of the SDL breaks the rotational symmetry and 
thus, on g eneral grounds, one expects the spiral core to 
move Il8j. This is indeed what is observed, albeit this 
motion is very slow, taking typically several thousand 
oscillations to move by one wavelength. After transients, 
the core moves ballistically at constant speed. In a fi- 
nite system, it eventually encounters a boundary which 
influences its motion. Figure El shows that the speed v 
increases continuously and monotonically from the onset 
of P2 bulk oscillations, and analysis of the data indicates 
a power law behavior v ~ {C — €2)^ , where C2 ^ 3.03 is 
the first period-doubling bifurcation and 7 = 1.5 ± 0.02. 



The angle a between the direction of motion and the at- 
tached SDL is constant (Fig. and gradually decreases 
from 180° to 90° with increasing C as shown in Fig. [21 

In the P4 and higher regimes both P2 and PI SDL may 
exist. We observed that both a and v are determined 
solely by the attached PI SDL and are unaffected by 
any P2 SDL that may be attached to the core, be it 
during a transient or in the regime where P2 SDL are 
spontaneously and continuously generated (4.306 < C < 
4.557). We conclude that if it exists at all, the effect 
of an attached P2 SDL is much weaker than that of the 
dominant PI line, and essentially undetectable. 

Next, we consider the turbulent regimes where PI 
SDL are spontaneously created (4.557 < C < 6.0). 
The continuous creation and annihilation of PI SDL 
and their dynamics strongly influences the motion of the 
core. More complicated configurations arise, which in- 
volve more than one PI SDL attached to the core, and 
the motion is no longer characterized by a simple angle a. 
Complex connection and reconnection events between PI 
SDL as well as between SDL and the core continuously 
occur. Close to the onset of this regime, the trajectory 
often changes its direction abruptly, even though long 
periods of (apparent) ballistic motion still persist. When 
a single PI SDL is attached to the core, the motion is 
ballistic with v > 0.0017 and a close to 90°. Yet, most 
of the time, the spiral core has three PI SDL attached to 
it (Fig. int) and short-lived configurations with a higher 
odd number of SDL exist as well. Even though a config- 
uration with three SDL is unstable for C < 4.557, it has 
a long lifetime and the combination of this persistence, 
the continuous creation of PI SDL, and the interaction 
with other nearby SDL (Fig.l^t), explains why this con- 
figuration dominates for 4.557 < C < 6.0. Interestingly, 
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Figure 3: Dynamics in the turbulent regime close to onset 
{C = 4.6, left) and far from onset {C = 5.8, right) for L = 512. 
(a): Snapshots of the Cx field, (b): Trajectories of the spiral 
core. The black arrow points to the position of the spiral 
core at the time when the snapshot of Acz was taken for 
C = 4.6 (shown in c). (c): Snapshot of Acz visualizing the 
line defects. Note that in the configurations shown three PI 
SDL are attached to the core. The white arrow points to the 
position of the spiral core. 



the direction of motion can change without any apparent 
variation in the configuration of the three attached PI 
SDL. For instance, during the change of direction high- 
lighted by the arrow in Fig. IHt), the core configuration 
shown in Fig. [Hh does not change significantly and the 
three SDL remain attached to the core. This indicates 
that part of the core motion is influenced by the dynam- 
ics of the background field in the turbulent regime and/or 
that subtle rearrangements occur deep inside the core. 

Deeper in the turbulent regime, the (apparent) ballistic 
flights become shorter and shorter and the core trajec- 
tory resembles Brownian motion (see Fig. IHb). In this 
regime, SDL are generated rapidly and homogeneously 
in the medium. This makes it possible to character- 
ize the core motion by a well-defined diffusion constant 



Figure 4: (color online) (a): Mean squared displacement 
(|Ar(t)|^) for different values of C in the turbulent regime, 
(b): The core diffusion coefficient Dy as a function of C. 



Dy. Figure ^ shows the mean-squared displacement 
(|Arv(t)n = (|rv(t) - rv(O)p) for different values of C. 
Here, (• • • ) denotes a time average over a long trajec- 
tory [JJI and Tv{t) is the position of the core at time t. 
The dependence of (|Arv(t)p) on t is clearly linear as 
expected for Brownian motion. The diffusion constant of 
the core motion determined from (|Arv(t)p) = ADyt de- 
creases with increasing C (see Fig.0lD), reflecting the fact 
that changes of direction become more frequent while the 
mean velocity remains roughly constant. 

We now turn to multispiral, spatially-disordered con- 
figurations which occur spontaneously in any large- 
enough domain. For oscillatory media described by the 
complex Ginzburg-Landau equation, it was shown re- 
cently that spatially-disordered frozen solutions do not 
exist I23. Even when the single spiral solution is stable, 
the weak but non-trivial effective interaction between spi- 
rals gives rise to ultra-slow core motion. These results 
should apply to the Rossler medium in the simple oscil- 
latory PI regime {C G [2.0,3.03]). A long run performed 
at C = 2.5 [isl showed that after a transient period a 
residual and barely detectable core motion subsists, with 
typical velocities decreasing sharply with increasing spi- 
ral domain size indicating a rapidly decaying interaction. 
This observation, together with the wide variation of spi- 
ral domain sizes (Fig. [5^), suggests that the PI Rossler 
medium is in the so-called "vortex glass" phase, and not 
the "liquid" phase where all spiral domains have roughly 
the same size. Thus, core motion exists also in this regime 
without SDL, albeit with velocities orders of magnitude 
smaller than those characteristic of core motion due to 
complex oscillations and can thus be neglected in the 
analysis of this latter case. 

In a spatial domain with periodic boundary conditions 
there must be an equal number of "positive" (clockwise) 
and "negative" (counter-clockwise) spirals. Under P2 lo- 
cal dynamics, spirals have a tendency to be grouped in 
pairs (mostly of opposite sign) connected by a PI SDL 
(Fig.Eb). Inevitably, the pattern rearranges since a con- 
nected pair drifts apart as shown in Fig. a single 
spiral pair in a square domain is connected by a straight 
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Figure 5: (color online) (a): Snapshot of Ac^ showing a 
typical multi-spiral configuration in the PI regime {C — 2.5, 
L — 1024). (b): same as (a) but in the P2 regime, where 
PI SDL connect spiral cores (C = 3.5, L = 512) (c): same 
as (b) but for a single pair of spirals connected by a PI SDL 
(L = 256). The direction of the cores' motion is highlighted 
by the white arrows, (d): Trajectories of spiral cores leading 
to the configuration shown in (b). Every 3000 time units the 
color changes starting with red or black depending on the 
orientation of the spiral. 

PI SDL. The two cores move with speed v at angles ±a 
relative to the SDL and drift away from each other with 
speed 2'L'cos(7r — a) which is positive as long as \a\ > 90°. 
The sign of a is determined by the sign of the spiral. The 
dependence of v and a on C is indistinguishable from the 
single-spiral case (Fig.|2|. 

In multispiral disordered P2 regimes our results show 
that spiral cores move nearly independently, stretching 
the PI SDL to which they are attached, until they 
meet another core leading to reconnections and cre- 
ation/annihilation of SDL. Some "frustration" arises dur- 
ing this process producing spiral cores connected by more 
than one SDL. The rearrangements occur mainly when 
two cores are close and the cascade of reconnection events 
can span a distance of several wavelengths. In a PI 
medium these cores would annihilate. On the contrary, 
here, because a > 90°, they repel each other after a new 
SDL connects them and the annihilation of spiral pairs 
is prevented (see Fig. Thus, multispiral configura- 

tions in P2 media are characterized by a novel type of 
spatiotemporal chaos which preserves the number of spi- 
ral cores through the complex dynamics of the PI SDL 
connecting them. 

Simulations show that this new "vortex liquid" state 



exists beyond the P2 regime, as long as a is significantly 
larger than 90°. For a close to 90°, annihilation of spiral 
pairs can occur, presumably due to the weaker repulsion 
effect induced by the PI SDL. Crossing the transition 
point to the turbulent regime at C = 4.557, the anni- 
hilation rate of spiral pairs drastically increases. This 
is caused by the continuous generation of new PI SDL 
which leads to permanent reconnections and effectively 
turns off the repulsion mechanism due to PI SDL. 

The new types of spiral core motion and turbulence 
described here should be observable in distributed me- 
dia with underlying complex dynamics. Studies of the 
more realistic Willamowski-Rossler reaction scheme 
support this statement. We also expect core motion to 
arise in excitable media supporting complex return to the 
rest state 22j. Experimentally, Belousov-Zhabotinsky 
reaction-diffusion systems are especially good candidates 
for studies of these phenomena since complex-periodic 
and chaotic regimes with SDL have already been ob- 
served in the laboratory |l5| . 
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